The Matter Bounce Curvaton Scenario 
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Massless scalar fields originating in a quantum vacuum state acquire a scale-invariant spectrum 
of fluctuations in a matter-dominated contracting universe. We show that these isocurvature fluctu- 
ations transfer to a scale-invariant spectrum of curvature fluctuations during a non-singular bounce 
phase. This provides a mechanism for enhancing the primordial adiabatic fluctuations and sup- 
pressing the ratio of tensor to scalar perturbations. Moreover, this mechanism leads to new sources 
of non-Gaussianity of curvature perturbations. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

A successful resolution of the hot Big Bang singular- 
ity may lead to a non-singular bounce cosmology. Non- 
singular bounces were proposed a long time ago 
They were studied in models motivated by approaches to 
quantum gravity such as modified gravity models @, Q , 
higher derivative gravity actions (see e.g. (IHsl), non- 
rclativistic gravitational actions 17]— [l Oi , brane world sce- 
narios [HEl, "Pre-Big-Bang" [ljand Ekpyrotic [3 
cosmology 1 , or loop quantum cosmology [18| . Bouncing 
cosmologies can also be obtained using arguments from 
superstring theory. For example, the string gas cosmol- 
ogy scenario may be embedded in a bouncing 
universe as realized in |2l| . Non-singular bounces may 
also be studied using effective field theory techniques by 
introducing matter fields violating certain energy condi- 
tions, for ex amp le non-conventional fluids I22L |23| , q uin- 
tom matter [2Ji28j , or ghost condensates [ly, [17], [29[ . A 
specific realization of a quintom bounce occurs in the Lee- 
Wick cosmology studied in (30 1. V arious original bounce 
models were reviewed in Ref. 31| . 

In the context of studies of bouncing cosmologies it 
has been realized that fluctuations which are generated 
as quantum vacuum perturbations and exit the Hubble 
radius during a matter-dominated contracting phase lead 
to a scale- invariant spectrum of cosmological fluctuations 
today [111, |33| (see also (3 for an earlier discussion). 
This yields an alternative to inflation for explaining the 
current observational data, which is dubbed as the matter 
bounce (see e.g. (35j for a recent review). However, in 
this scenario there are still unclear issues. For example, 
since scalar and tensor modes grow at the same rate in 
a matter-dominated phase of contraction, the tensor-to- 
scalar ratio is typically too large to be consistent with 
current observational data. Also, the possible generation 
of entropy during the non-singular bounce phase needs 



to be studied in detail. In this paper, we will focus on 
the first question, leaving the second for a followup paper 

M- 

To date, entropy fluctuations have not been considered 
in matter bounce models. However, realistic matter the- 
ories will contain many entropy modes. In the context of 
inflationary cosmology, it has been realized that entropy 
fluctuations can lead to an additional source of curvature 
fluctuations, a source which may in fact dominate 
l4lj . This is now know as the "curvaton mechanism" for 
generating fluctuations. Entropy fluctuations also play 
an important role in the Ekpyrotic scenario 
13- I n that scenario, the primordial adiabatic fluctua- 
tions in the contracting phase acquire a steep blue spec- 
trum, whereas the entropy modes induced by light scalar 
fields will be scale-invariant. 

In this paper, we will consider the fluctuations induced 
by a light spectator field x an d see that they acquire 
a scale-invariant spectrum of fluctuations in a matter- 
dominated phase of contraction if they originate in their 
vacuum state on sub-Hubble scales. We will show that 
these entropy fluctuations induce a scale-invariant spec- 
trum of curvature perturbations after the bounce. The 
transfer of entropy to adiabatic perturbations occurs dur- 
ing the bounce phase. We can easily find parameter val- 
ues for which these secondary curvature perturbations 
are larger in magnitude than the primordial ones. This 
will lead to a suppression of the tensor to scalar ratio. 

As a simple example, we consider a two-field model 
with potential 



V(cf>, X ) = \m 2 <P + i 9 W 



(1) 



where in our model </> is the background scalar field in the 
nonsingular bounce cosmology, and x is a second scalar 
matter field, the so-called entropy field 2 . 



1 Some versions are non-singular, see e.g. |15H in the case of the 2 The non-singular bounce is induced by adding another matter 
Pre-Big-Bang scenario, or the "New Ekpyrotic scenario" [l6l , ll7ll . field violating the null energy condition to the system, or by 
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In the context of a matter bounce background, we 
study the cosmological evolution for both background 
and perturbations. We study the specific model given 
by the potential dJ). We find that the entropy field x 
evolves as a tracking solution in the contracting phase 
due to the coupling to the background scalar </>, and may 
become dominant during a short deflationary phase be- 
fore the bounce. The fluctuations of the entropy field 
are nearly scale-invariant provided the coupling param- 
eter g is small enough. A controllable amplification on 
the entropy modes due to the kinetic term of the entropy 
field can be obtained in the bouncing phase, and these 
entropy modes will be converted to curvature perturba- 
tions during the bounce. Therefore, we can obtain an 
enhancement of the primordial curvature spectrum and 
hence a suppression of the tensor-to-scalar ratio. This is 
the main message of this work. Another interesting fact 
is that the conversion of entropy fluctuations to curva- 
ture perturbation provides a sizable non-Gaussianity for 
curvature perturbations for which the sign tends to be 
negative. 

The paper is organized as follows. In section II we 
provide the basic equations of the theory. In section III 
we study the generation of isocurvature fluctuations in 
the contracting phase. In section IV we analyze the con- 
version of the isocurvature perturbations into curvature 
fluctuations during the bounce and show that a reason- 
able amplitude of the primordial power spectrum can be 
obtained. Section V presents a summary and discussions 
of related works. 

We use the convention m p i = 1/yG in this paper. 



II. SETUP AND BASIC EQUATIONS 

To be specific, we will be considering the background 
space-time which emerges from a quintom bounce (see for 
example [30| ) scenario in which regular matter is mod- 
elled by a scalar field <fi with mass m and the non-singular 
bounce is induced by adding a ghost scalar field ip g with 
a much larger mass M to the system. Fig. [T]is a sketch of 
the evolution of scales in our matter bounce model. The 
vertical axis represents conformal time, the horizontal 
axis comoving distance. Initially, both the regular matter 
field and the ghost field arc oscillating, with amplitudes 
which are increasing since the universe is contracting. 
During this period, the Hubble radius is decreasing lin- 



considcring the matter Lagrangian containing the regular fields 
<j> and x m the presence of a gravitational action leading to a 
bounce. Note that in a successful bouncing cosmology the usual 
energy conditions need to be recovered after the bounce so that 
the universe is able to evolve smoothly into the radiation and 
matter dominated periods. This feature exactly coincides with 
what occurs in the quintom dark energy model 14511 ■ in which the 
dark energy equation of state crosses the cosmological constant 
boundary. 



early and H < 0. We assume that the energy density is 
dominated by the regular field <p- Once the amplitude of 
the regular field (j> reaches the value m p i, the field will 
cease to oscillate. This happens at conformal time r m . 
A short deflationary "slow climb" phase will begin dur- 
ing which the Hubble radius \H | _1 is nearly constant. In 
this period the ghost field continues to oscillate. Hence, 
the (negative) energy density of the ghost field catches 
up to the (positive) energy density of the field <f> and a 
short non-singular bounce period will result during which 
H > 0. The conformal time is the transition between 
the deflationary phase and the bounce period. After the 
bouncing phase there is a short period of inflation lasting 
until tri after which the universe enters the expanding 
phase with a normal thermal history. We denote the 
bounce time by tb- 
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FIG. 1: A sketch of the evolution of scales in a nonsingular 
bouncing universe. The horizontal axis is comoving spatial 
coordinate, the vertical axis is conformal time. Plotted are 
the Hubble radius \H | _1 and the wavelength A of a fluctuation 
with comoving wavenumber k. 

In Fig. Q] we plot the evolution of the length corre- 
sponding to a fixed comoving scale. This scale is the 
wavelength of the fluctuation mode k (k standing for the 
comoving wavenumber) which we want to follow. The 
wavelength begins in the contracting phase inside the 
Hubble radius, exits the Hubble radius during this phase 
and re-enters the Hubble radius during the normal ex- 
panding phase. 

It is important that scales of interest exit the Hub- 
ble radius in the matter-dominated phase of contraction. 
This is required to obtain a scale-invariant spectrum of 
cosmological perturbations. Important is also that there 
is a phase during which the energy density in the entropy 
field grows relative to the energy density in the field <f>. 
This condition is easily satisfied in the quintom bounce 
scenario. It is required in order that the conversion of 
the entropy fluctuations to curvature perturbations is ef- 
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ficient 3 . The efficiency of our transfer mechanism will 
depend on whether the initial energy density of the en- 
tropy field x is larger than the (absolute value of) the 
energy density in the ghost field ip g . We shall see that 
our mechanism is more efficient if \ dominates 4 . 

The key goal of our paper is to understand how entropy 
perturbations arc generated in the contracting phase, 
then pass through the bounce, and finally convert into 
curvature perturbations after the bounce. In the follow- 
ing, we will study cosmological perturbation theory in the 
nonsingular bouncing cosmological background described 
above. 

We shall work in longitudinal gauge in which the lin- 
earized scalar metric fluctuations appear in the metric in 
the following way: 



ds 2 = (1 + 2$)d< 2 - a 2 (i)(l - 2$>)da 



(2) 



where t is cosmic time and x l are the comoving spatial 
coordinates. The scalar metric fluctuations are charac- 
terized by two functions $ and \t which depend both on 
space and time. We take matter to consist of a set of 
scalar fields fa. If the gravitational action is the usual 
one, then in the absence of linearized anisotropic stress 5 
the off-diagonal components of the perturbative Einstein 
equations imply ^ = $ (see j47| for a comprehensive dis- 
cussion of the theory of cosmological perturbations). By 
expanding the Einstein and matter equations to first or- 
der, we obtain the following perturbation equations (48l — 



3H5fa 



V 2 



= 40i$ - 2V,i$ , 
V 2 

-3iJ$ + (— -3ff 2 )$ 



(3) 



= iirG^lMi - </> 2 $ + VJfa] , (4) 

i 

<& + H<f> = inG^faSfa , (5) 



where H = a /a is the Hubble parameter, and Vi denotes 
the derivative of the scalar field potential with respect to 



3 This is similar to the requirement in inflationary cosmology for 
the curvaton mechanism to be efficient. 

4 It would be of great interest to study to what extent the conver- 
sion of entropy fluctuations to curvature perturbations is efficient 
in other bounce models such as the ghost condensate scenario of 
|29| which is (unlike the quintom bounce 46]) free of ghosts at 
the perturbative level, and stable against the addition of radia- 
tion and anisotropic stress. 

5 This will be the case if matter is described by a set of scalar 
fields. 



We can recast the above equations in terms of the 
Sasaki-Mukhanov variables [fLU, which are defined as 



Qi = Sfa + , 



(6) 



and in terms of which the equations of motion are given 
by 

Qi + 3HQi - ^-Qi + X^y - ^-{^fah)]Qi = • 

3 

(7) 

To combine the above equations, we usually define the 
quantity C which is the curvature perturbation on the 
uniform density slice, and which is given by 



C = H 



J2i <PiQi 



(8) 



This quantity is conserved on super-Hubble scales in an 
expanding universe if there are only adiabatic fluctua- 
tions. However, the presence of entropy fluctuations on 
large scales will lead to a growth of ( which corresponds 
to the seeding of an adiabatic fluctuation mode by the 
entropy mode. 

At linear order, the general equation for the time 
derivative of £ is (see Appendix) 



H 4^ 6 4 H J 



(9) 



On large scales, the first term of the r.h.s of Eq. ([9]) is 
negligible. In a two field model, this equation reduces to 
the following expression for the time derivative of Q on 
large scales: 



(l + w)C = AttG- 



+ X 2 <Q4> 



3H 



X dt 02 + x 



77) . (io) 



where w = p/p is the equation of state, with p and p being 
pressure and energy density respectively. Note, however, 
that in our case we have three fields, the regular matter 
field (j), the isocurvature field x an d the ghost matter 
field tpg which is responsible for yielding the non-singular 
bounce. 

In the following we will study how in our non-singular 
bouncing cosmology the presence of primordial isocur- 
vature fluctuations induces via ^ a contribution to the 
curvature perturbations. We will consider two cases: the 
first where the isocurvature field x initially contributes 
more to the density than the ghost field, the second case 
where x is sub-dominant at the background level to the 
ghost field. 



III. ISOCURVATURE PERTURBATIONS IN 
THE CONTRACTING PHASE 

In this section we show that the isocurvature mode ac- 
quires a scale-invariant spectrum of fluctuations on scales 
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which exit the Hubble radius in the matter-dominated 
phase of contraction, provided the initial conditions are 
taken to be vacuum. We will consider the model of ([1}, 
where the entropy mode corresponds to the scalar field \ 
which is taken to be massless. Note that since the ghost 
field ip g has a large mass, its spectrum will not be scale- 
invariant. It will be highly suppressed in the infrared. 
Hence, ip g cannot play the role of the entropy field for 
our purpose. 

We consider a matter bounce background evolution of 
a FRW universe in which the matter-dominated phase 
of contraction is realized by the oscillations of a mas- 
sive scalar field (f> around it vacuum state (see Rcfs. 
[Ill, HH, HH, Hil for other realizations). Hence, the Hubble 
parameter can be approximately expressed by 



acquire a scale-invariant spectrum in the matter bounce 
scenario is based on assuming that the canonically nor- 
malized curvature fluctuation variable starts out in its 
vacuum state early in the contracting phase. It is hence 
reasonable to assume that also the fluctuations of x start 
out in their vacuum state on scales smaller than the 
Hubble radius early in the contracting phase. We fo- 
cus on scale which exit the Hubble radius in the matter- 
dominated phase of contraction and will show that due 
to the growth of fluctuations on super-Hubble scales a 
flat spectrum of isocurvature perturbations results 6 . 

We are interested in the scalar perturbation 6x, where 
we write 



X (t,x) = x(t) + Sx(t,x) 



(16) 



(H(t)) 



2 

M 



(11) 



where the angular brackets stand for averaging over time. 
The solution for the scalar field 4> can be asymptotically 
expressed as 



cj)(t) ~ (f>(t) smmt , 
with a time dependent amplitude, 



37rm|i| 



(12) 



(13) 



which yields an equation of state which has vanishing 
pressure after averaging over an oscillation period of the 
field. 

Since the entropy field x is much lighter than the back- 
ground scalar <f> and we assume that it does not have a 
large initial homogeneous mode, it has nothing to do with 
the background evolution before the bounce phase. Its 
evolution is determined by the equation of motion 



X 



3H X + ^9 2 4> 2 X - , 



(14) 



where we have averaged the coherent oscillations of the 
background scalar in the mass term and inserted a factor 
1/2 instead (parametric resonance instabilities resulting 
as a consequence of the oscillatory nature of <j> (55l - [58j - 
see [59j for a recent review - will be studied in a subse- 
quent paper). This equation yields the following solution 
for x- 



X(t) oc t 



-1+'- 



(15) 



which evolves at the same rate as the amplitude of the 
background scalar if gm p i <C m. In this case it is conve- 
nient to define a parameter d± = & (where the subscripts 
"i" indicate the initial value of the corresponding fields), 
and then in the matter-dominated contracting phase we 
have x — dx4>. 

Let us now turn to the fluctuations of x- Let us re- 
call [H, HH that the finding that curvature fluctuations 



It is important to focus on the canonically normalized 
perturbation variable v = a5x in terms of which quan- 
tum vacuum initial conditions can be imposed. We work 
on the spatially-flat slice where the metric fluctuation 
$ vanishes. On this slice, the scalar field perturbations 
satisfy the Klein-Gordon equation 61 



< + (r + am H 

* a 







(17) 



where we have defined the effective mass term of the en- 



V xx and made use of the conformal 



tropy field m x 
time drj = adt. A prime indicates a derivative with re- 
spect to rj. Comparing the above equation of motion with 
the basic perturbation equation ([7]), one may notice that 
the effective mass originating from the kinetic term of x 
and the coupling between the isocurvaure and adiabatic 
perturbations has been neglected. This is justified since 
it is subdominant in the matter-dominated contracting 
phase. 

Initially, the k 2 term dominates Eq. (fTT)) . and so we 
can neglect the mass term and gravitational term. Thus, 
the dynamics corresponds to a free scalar propagating in 
a flat spacetime, and the initial condition takes the form 
of the Bunch-Davics vacuum: 



-ikr 



2k 



(18) 



During the phase of background contraction, the quan- 
tum fluctuations exit the Hubble radius and become clas- 
sical perturbations (see [62|, [63| for an analysis of the 
classicalization). For a oc t 2 and \m 2 x \ < H 2 , Eq. ([T7]) 
becomes 



(k 2 



rH - 



(19) 



6 As discussed in [(Still , the spectrum has a break at the scale which 
exits the Hubble radius at the time of the transition between the 
matter-dominated phase of contraction and the following phase 
(most likely a radiation-dominated phase). On smaller scales the 
spectrum will be blue. 
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The second term inside the parentheses on the left hand 
side of this equation is the term which leads to the squeez- 
ing of the fluctuations, and the coefficient and time de- 
pendence is exactly that required to transform a blue vac- 
uum spectrum into a scale- invariant spectrum [HI, |33| . 

We can see this explicitly as follows: Making use of the 
vacuum initial condition, we obtain an exact solution to 
(HI: 



-ikr 



<■/,■ 



(i 



kr' 



(20) 



Note that, in contracting phase these perturbations are 
able to grow outside the Hubble radius. Their ampli- 
tudes will keep increasing until the universe enters the 
contracting phase ends. Thus, from (|20[) it follows that 
on scales larger than the Hubble radius, the amplitude of 
the spectrum of Sx in the contracting phase is given by 



\H(t)\ 

47T 



(21) 



which is a function of cosmic time. From Eq. (|17j) . we 
are able to calculate the spectral tilt n x of the primordial 
perturbations 



d\nP x 
din k 



2m 2 x 
3#2 



(22) 



where m x is the effective mass of the x held induced by 
the coupling to <j>. In the specific model given by ([1]), we 
apply the amplitude of the background scalar <f> and the 
ratio d, and then obtain the spectral tilt 



2 2 

9 m P i 
27rm 2 



(23) 



From this result, we find that in order to ensure that the 
isocurvature perturbation are nearly scale-invariant, one 
has to require |<77n p ;| <C \phim. 

One can also estimate the amplitude of tensor pertur- 
bations h. It is well known that the tensor fluctuation 
amplitude satisfies the same equation of motion as a free 
scalar field in the given cosmological background. Hence, 
following the analysis done above for the fluctuations in 
the scalar field Xi we obtain a scale-invariant spectrum 
of tensor modes (again assuming that they start out in 
their vacuum state), and an amplitude of the spectrum 
given by 



H/m p i . 



(24) 



The tensor perturbation amplitude has the same order 
of magnitude as that of the adiabatic scalar metric per- 
turbations. Thus (in terms of the usual ratio r = Pt/ ' Pq 
of the tensor to scalar metric power spectrum) , the mat- 
ter bounce usually leads to a dangerously large value of 
r. From current observational CMB data [64| . this ratio 
is required to be less than 0.2. Consequently, in order 
that a matter bounce scenario be consistent with current 



observations on the value of r, it is necessary that there 
be a mechanism to amplify the size of the scalar met- 
ric perturbations. In the following we will show that the 
entropy mode studied in this paper can provide an effi- 
cient mechanism to enhance the spectrum of curvature 
fluctuations. 

The basic idea is as follows: We expect that the en- 
tropy field will rapidly grow during the final instances of 
the contracting phase. In this case, isocurvature pertur- 
bations can be converted into curvature perturbations as 
described by ([TU| 7 . After the bounce, the evolution of 
the universe will enter the normal thermal history. In 
contrast to the bounce scenario without entropy modes, 
in the current model the energy scale for the bounce is 
able to be lowered and so the amplitude of tensor spec- 
trum can be suppressed. Therefore, our model can solve 
the problem of large tensor-to-scalar ratio which gener- 
ally exists in nonsingular bounce cosmologies (30 . l6lj . 
Moreover, if the energy scale can be lowered to the TeV 
level, our scenario may lead to a window to test a bounce 
in particle physics. 



IV. EVOLUTION OF ISOCURVATURE 
PERTURBATIONS THROUGH THE BOUNCE 

As the universe contracts in the matter-dominated 
phase, the Hubble parameter increases. When it catches 
up to the mass m of the background field (f>, the oscilla- 
tions of this field will freeze out. After freeze-out, <fi will 
slowly climb up the potential. As long as the total en- 
ergy density is dominated by (f>, the equation of state will 
decrease from w = to w ~ - 1. Therefore, the universe 
will enter a contracting deflationary period and then ap- 
proach the bouncing point. During this phase, both the 
isocurvature field x an d the ghost field tp g will increase 
exponentially in amplitude. Thus, after a brief deflation- 
ary phase either x or "4>g will come to dominate. If the 
energy density of ip g is larger, then the bounce point is 
reached when the energy densities of <p and ip g become 
equal. We will denote this case as "Case 2". In the 
other case ( "Case 1" ) there will be another brief phase of 
contraction driven by Xj followed by another brief defla- 
tionary phase and a bounce when p x = |p-i/>J. Note that 
the deflationary phases can be understood as the time 
reversal of a period of slow-roll inflation. 



A. Deflationary period 

We first study the evolution of the entropy scalar x 
in the deflationary period. From studies of chaotic type 



7 As we discuss in our companion paper [3611 . they can also be 
paramctrically amplified when they pass through the bouncing 
phase. 
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inflation models, it is well known that the background 
field </> freezes out when its amplitude <j) becomes of the 
order of the Planck mass, more specifically when 



in Case 1 and 



2m r , 



This leads to a nearly constant Hubble parameter 
rr - — 



(25) 



(26) 



and an almost exponentially decreasing scale factor for 
the universe. In this period, the scalars \ an d ipg will 
grow exponentially and the density in one of them will 
soon catch up to the density of background field. The 
background equation of motion (|14[) yields 



X 



Xm exp 



4m 



(t - t m ) 



(27) 



which evolves proportional to a~ 3 . Here, the time t m is 
the end of the matter-dominated epoch of contraction, 
and the subscripts "m" on the fields <j> and x indicate the 
values of those field at the time t m . 
We will write 



= d 2 4> m 



(28) 



where the values of the coefficients d\ and d 2 (both much 
smaller than 1) depend on the initial conditions. We will 
denote the larger of the two numbers d± and di by d. 
The deflationary phase will end when the energy of one 
of the fields x an d ip g catches up to that in <p. This will 
occur at the moment t^ given by 



td — t, 



7T j 7T 

4m 4di 



in Case 1 and 



7r m 



(29) 



(30) 



in Case 2. The difference between the two cases is due to 
the fact that ip g is oscillating in the deflationary phase, 
but x is n °t- Hence, the expressions for the energy den- 
sity of the two fields take a different form. 

Let us now turn to the evolution of the entropy fluc- 
tuations which is given by Eq. (|T7|) . In a de Sitter phase 
of contraction the dominant solution is 



(34) 



in Case 2. 

Note, however, that this amplification does not lead to 
a suppression of the tensor to scalar ratio since the ten- 
sor modes at the same rate during the deflationary phase. 
It is only due to matter couplings that a difference be- 
tween scalar and tensor modes will result. Below, we will 
see that such an amplification mechanism affecting only 
the scalar modes results in the post-deflationary phase, 
a phase driven by the kinetics of the fields. 

We do need to worry about the absolute magnitude of 
the fluctuations. There is the danger that the amplifi- 
cation in the deflationary phase cause the perturbations 
become nonlinear. Then, we would lose perturbative con- 
trol of the analysis. Presumably, this would also imply 
that space fragments into a gas of black holes rather than 
undergoing a smooth bounce. However, from (|33|) and 
from the fact that the fluctuations start out as quantum 
vacuum perturbations, it appears that this constraint on 
the parameters is a fairly mild one. 



B. Passing through the bouncing phase 

In this subsection we will give a rough analysis of how 
the fluctuations in x pass through the bouncing phase. 
We will discover that the coupling between the matter 
fields due to the kinetic terms results in a growth of the 
entropy fluctuations without corresponding growth of the 
tensor modes. We call this process kinetic amplification. 

After the period of deflation, the scalar x continues to 
roll fast, and may even dominate over the background 
scalar <f> (in Case 1 - see the beginning of this section 
for the definition of what the Cases are). In this case, 
the universe evolves with a background equation of state 
which varies from w ~ — 1 during deflation to w ~ 1 
during x domination. In this case, the dominant term in 
the background Friedmann equation is the kinetic term 
of the scalar \i an d leads to a nearly conserved quantity 



X/H ~ Xd/Hd ^ Xd/Hr, 



(35) 



during this short period. In addition, from Eq. (|27l) wc 
see that 



Xd 



4m 4m ~ 

Xd ^ <Pr. 



(36) 



which corresponds to 



o~Xk 



(31) 



(32) 



Hence, it follows that the power spectrum of 5x gets en- 
hanced during the deflationary phase by a factor 



C 



( — ) 2 



(33) 



As discussed earlier, eventually the energy density in 
the ghost field ip g catches up and the bounce phase will 
be triggered. Thus, the universe will exit the contraction 
phase and the equation of state of the universe will cross 
— 1 and fall to negative infinity. In order to study the 
evolution of perturbations through the bouncing phase 
analytically, it is convenient to model the evolution of 
the Hubble parameter in the following form: 



H = at 



(37) 
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with some positive coefficient a which has dimensions of 
k 2 and whose magnitude is set by the microphysics of the 



bounce. For example, a 



in the model of Lee- Wick 



Cosmology [30|, while in Loop Quantum Cosmology a is 
determined by a regularized Planck mass[65l fr36l|. Note 
that in the above paramctrization we have chosen the 
time of the bounce to be zero. 

To study the perturbation of the \ field, it is necessary 
to go back to the complete form of the perturbation equa- 
tions given by ([7]). We will focus on the evolution of Q x . 
From ([7]) we see that there is a mixing with Q^. However, 
since the frequency of 8<f> is dominated by its mass to, 8<p 
keeps oscillating even outside the Hubble radius. Hence, 
its amplitude is much smaller than that of 8\ an d there- 
fore we neglect its contribution in the equation of motion 
for 8x- Thus, Q x effectively decouples from Q^. Focus- 
ing now on the mass term in the equation of motion for 
Q x , we notice that since x evolves as a stiff fluid with an 
effective equation of state w = 1, the contribution of its 
kinetic term in the perturbation equation must be taken 
into account. This contribution is given by the following 
term: 



a 3 m 2 pl \ H 



8tt /a 3 x 2 



727T02 



(a + 3arr) , (38) 



where we have used the relations (|35|) . (|36|) and (f37|) . 

Consequently, we obtain the following equation of mo- 
tion for the Fourier modes of the gauge invariant pertur- 
bation variable Q x , 



« 2 



Q x + 3HQ X + '—Q x + m 2 x cS Q x ~ . (39) 



In the bouncing phase the effective mass term is given by 



l x,eff 



212 
9 <Pm ~ 

4g 2 m 2 



72naq 



Pi 



3tt 3 



"V 
96a 



-(l + 3crf 2 ) 



3ai 2 ) 



(40) 



Thus, we find that in the bouncing phase the contribution 
of the kinetic term to the perturbations of the x field 
yields a tachyonic mass term. This results in an increase 
in the magnitude of the fluctuations, a process which we 
call kinetic amplification. Note that this amplification 
does not affect the tensor modes. 

In analogy with the theory of preheating (see e.g. [59| 
for a recent review), we suggest that the investigation of 
the perturbation Q x can be simplified by redefining the 
fluctuation variable 



,3/2 



Q 



Xk 



(41) 



Then, we obtain a simpler form of the equation of motion 
for the fluctuations. Instead of Eq. (|3"9")) we now have 



SXi. 



isx k = 



(42) 



where 



J k = — 



l X ,eff 



3 • 



(43) 



is the frequency for the k'th mode of the isocurvature 
perturbations. 
Consequently, if 



a > g 2 m 2 pl , 
there exists a period — i& < t < tf, with 
h ^ \H m \/a = Am/San 



(44) 



(45) 



near the bounce point during which the adiabaticity con- 
dition is violated and the dispersion relation of the x per- 
turbation is unstable in the infrared limit. In this case 
the evolution of the x perturbation can be approximately 
expressed as [63] 



SX oc exp[| J dtuj(t)\] 



(46) 



in which the integral runs from — tb to tb if a 3> g 2 m'^ )l . 

As a consequence, we find that the amplitude of the 
entropy fluctuations obtains a tachyonic amplification 
around the bounce point. Since this instability arises 
from the kinetic term of the entropy field, we call this 
process kinetic amplification. One can describe the am- 
plification of the infrared modes of isocurvature pertur- 
bations in this resonance stage as follows: 



Sxb = FSxd , 
with an amplification factor 

/tb 
dtw(t)\] 
-tb 



(47) 



T 



v /^Ts0+4 r sinh- 1 (i£2 



(48) 



In the above calculation, we have introduced a variable 

2 

y = — which characterizes how fast the bouncing phase 
occurs. 

To understand the exponent appearing in the ampli- 
fication factor J-, we plot its evolution as a function of 
the variable y in Fig. [2] One can see that in the limit 
of a fast bounce y — > (which corresponds to a > m 2 ) 
the contribution of kinetic amplification is quite limited 
since the duration of the bouncing phase is very short. 
In the case of the Lee- Wick cosmology a ~ to 2 and thus 
gives y ~ 1, and correspondingly we obtain a value of the 
amplification factor J- ~ 21.3 which is quite sizable. 

One interesting issue which should be noticed that, if 
we choose the value of the coefficient a to be between 
y 2m pi ano - m2 ' the- amplitude of the isocurvature pertur- 
bation will be distinctly amplified. This case corresponds 
to a slow bouncing phase. This amplification becomes 
more efficient if we take the coupling parameter g 2 to be 
slightly negative which provides a flat tachyonic potential 
for the entropy field. 
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FIG. 2: Plot of the exponent in the amplification factor XnT 
as a function of the variable y in the bouncing phase. The 

2 

permitted range of y runs from to , m a , and we specifically 

9 ~ m v l 

choose g 2 m,p l = m 2 /16. 



C. The x power spectrum in the expanding phase 

After the bounce, the universe will eventually re- 
cmcrge into a matter-dominated phase, this time a phase 
of cosmological expansion. Since the dominant mode 
of both scalar and tensor fluctuations on super-Hubble 
scales is constant, no further amplification of the entropy 
fluctuations results, nor do scalar and tensor perturba- 
tions evolve asymmetrically. 

Now let us take into account both of the amplification 
effects on the x fluctuations discussed in the two previous 
subsections. Combining the relations ([55]) and (|4"T)) and 
applying the relation (|21j) at the beginning of the period 
of deflation at the time t m , we eventually obtain the fol- 
lowing amplitude of the isocurvature perturbations after 
the bounce 

S X , - , (49) 

on super-Hubble scales. 



V. GENERATION OF METRIC 
PERTURBATIONS 

In the above section we have studied the dynamics of 
isocurvature perturbations in a particular implementa- 
tion of the matter bounce scenario. Next, we will study 
the induced curvature fluctuations, the "matter bounce 
curvaton" mechanism. The transfer takes place in the 
deflationary phase and in the bouncing phase when the 
equation of state of the background is changing because 
of the contribution of the entropy field. 



A. Linear fluctuations 

We first study the linear curvature fluctuations induced 
by the isocurvature fluctuation. To do this, we return to 
the general formula ^ which describes the transfer of 
isocurvature to curvature fluctuations. Making use of 
the fact that 

m pi i 

and inserting the definition of Qi we find that on large 
scales ([9]) reduces to 

Since the spectrum of the ghost field ip g has a blue 
spectrum due to the large mass of the ghost, the ghost 
field contribution to ([51"]) can be neglected, and only two 
terms remain in the sum, the one due to 4> and the one 
due to the entropy field. Since in our realization of the 
bounce scenario the regular matter field <f> is massive, its 
contribution to ([5Tj) is also negligible. A second way to 
see that the term proportional to S(f> gives a negligible 
contribution to (|51|) is to realize that until the entropy 
field x begins to dominate the background, the factor 
4> 2 /H is approximately constant and thus the coefficient 
of the 8(j) term vanishes. Thus, (|51[) can be approximated 
by 

C - ■ OB) 

X ti 

The right hand side of this expression vanishes at the 
beginning of the deflationary phase and gradually grows 
until the bounce phase is reached. We can approximately 
integrate the equation to obtain 

AC ~ ^nm-^xSx- (53) 

Let us first evaluate the magnitude of the induced cur- 
vature fluctuations in Case 1, when p x > \pi/) g \- In this 
case the contribution of the ghost condensate field to H 
at the time when the deflationary period ends can be 
neglected and 

H ~ -47rm- 2 (0 2 +x 2 ). (54) 

When the x starts to dominate the energy density, 
the magnitudes of x an d <j> are comparable, and hence 

AC ^ -5 X . (55) 

X 

Making use of the background solution for the scalar x 
and the amplification of the amplitude of the isocurva- 
ture perturbation 5xf before and during the bounce, we 
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finally obtain the following estimate for the induced cur- 
vature perturbation 



AC 



6v37rm p / 



(56) 



where we have used relation (|35[) which is approximately 
valid during the deflationary and bounce periods. 

Next, let us turn to Case 2. Now, H is dominated by 
ip g instead of by \ m t ne bounce region. The numera- 
tor in ([53]) remains the same as for Case 1. Hence, the 
overall amplitude of the induced curvature fluctuations 
is suppressed by a factor S, where 



S 



(f) 2 , 



(57) 



compared to what happens in Case 1. 

The tensor modes are also amplified during the de- 
flationary phase, as mentioned earlier. Thus, they also 
obtain the amplification factor C. Their overall ampli- 
tude is hence the same as that of the induced curvature 
fluctuations, except for the factor T due to kinetic am- 
plification. Hence, the predicted tensor to scalar ratio r 
in Case 1 is 



Our analytical calculations involve some approxima- 
tions. Specifically, in the deflationary phase the Hubble 
parameter is not exactly a constant and the background 
scalar <p is still rolling slowly, but in our analytical analy- 
sis we have taken their values at the freeze-out time. As 
a consequence, it is important to confirm the results by 
a numerical integration of the full equations. 

We first numerically study the evolution of two scalars. 
As demonstrated in Fig. [3l in the era of contraction, the 
heavy field <p keeps oscillating with an increasing ampli- 
tude and dominates the background evolution, while the 
light field x tracks the evolution of <p. When the am- 
plitude of 4> grows up to of Planck mass, x rolls rapidly 
and finally catch up to cp near by the bounce. This nu- 
merical result is consistent with the analytic estimates 
appeared in Eqs. (|15[) . (|27[) . However, in our numerical 
calculations the background evolution ceases around the 
moment of x domination which is before the big bounce. 
Since in the current paper we wish to discuss the dy- 
namics of isocurvature perturbations in a general case, 
we did not concentrate on a detailed mechanism leading 
to a bouncing phase but take a parametrization of the 
Hubble parameter H — at instead. 



(58) 



Its magnitude is determined by two parameters, namely 
the mass of the heavy field m and the slope of the bounc- 
ing phase a. 

According to the latest observations from the Wilkin- 
son Microwave Anisotropy Probe 7- year data (WMAP7) 
[64j| . the amplitude of curvature perturbation is C — 
4.92 x 10~ 5 and the spectral tilt is \n x \ < 0.03. Thus, 
the observed amplitude of the curvature perturbations 
can be obtained in our scenario if we require 



9.02 x 10~ 4 J^TOpz 



(59) 



For completeness, let us mention (as was already used 
above in (1581) ) that in our model the amplitude of tensor 
perturbation is about h ~ H/m p . Specifically, following 
the calculations in Ref. [3(|, one obtains the following 
primordial power spectrum for the tensor perturbations 
at the end of the matter dominated phase of contraction 



P T (k) 



32k 3 



\h k \ 2 



16H 2 
9Trm 2 pl 



(60) 



During the deflationary phase, there is an enhancement 
by a factor of C. The tensor-to-scalar ratio is given by 




FIG. 3: Evolution of the background fields cj> and \ before 
the bounce as a function of cosmic time (horizontal axis). 
The background fields are plotted in dimensionless units by 
normalizing by m p ;. Similarly, the time axis is displayed in 
units of m~; . The parameters m and g appearing in the 
potential were chosen to be m = 10 _5 m p ;, g = 2.5 x 10~ 6 , 
and we assume initially d = 10~ 4 . The initial conditions for <j> 



were taken to be fa = l.I9x 10 5 m p i, fa = — 1.24x10 9 m p i 2 . 



35J" 



(61) 



This ratio can be greatly suppressed by a large value of 
the factor JF . In the specific model of Lee- Wick cosmol- 
ogy, a ~ m 2 so that r ~ 0.08. With the WMAP7 data 
combined with BAO and SN, the current limit on the 
tensor-to-scalar ratio is given by r < 0.2. 



Based on the background solution obtained above, we 
then solve the equations of motion <[7|) numerically. We 
show the evolution of the gauge-invariant scalar pertur- 
bation variables and Q x in Fig. 2J In the numer- 
ical calculations, we consider these perturbation modes 
to originate from vacuum fluctuations inside the Hubble 
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radius 8 . One can see from Fig. 2] that the perturbation 
of the heavy field (the adiabatic mode after the bounce) 
catches up to the amplitude of the entropy field perturba- 
tion. This demonstrates the effectiveness of the "matter 
bounce curvaton" mechanism discussed here. Note that 
Qcf, keeps oscillating through the whole evolution. The 
amplitude of the fluctuation Q x grows before the bounce, 
however, with a different slope than the growth of Q^. 
This difference in slope is due to the transfer of isocur- 
vature to curvature fluctuation. Overall, the numerical 
results are consistent with the analytic calculations. 




FIG. 4: Evolution of the gauge-invariant scalar perturbations 
Q</> and Q x for super-Hubble mode k — 2 x 1(F 6 in our model. 
The horizontal axis is cosmic time. The scalar perturbations 
are plotted in dimensionless units by normalizing by m p i. In 
the parametrization function a = 10 _e m?j is chosen. The 
background parameters are the same as in Fig. 



B. Local non-Gaussian fluctuations 

Next we take a first look at an important byproduct of 
the "matter bounce curvaton" scenario, namely a contri- 
bution to the induced non-Gaussian distribution of met- 
ric fluctuations which will be of local form. 

There are two important periods for generating a large 
non-Gaussianity in a bouncing universe. The first is the 
evolution of metric perturbations on super-Hubble scales 
in the contracting phase which at leading non-linear order 
generates non-Gaussianities as studied in Refs. [zHllBI- 
This generates an effective value of the non-linearity pa- 
rameter f n i which is of order one, the largeness being 



due to the fact that what corresponds to the inflation- 
ary slow-roll parameter is now a constant of order one. 
A distinctive feature of this source of non-linearity is a 
particular shape of the bispectrum. The shape is dif- 
ferent than the one which appears in simple inflationary 
models since in the contracting phase £ is not constant 
and hence terms depending on £ which vanish in infla- 
tionary cosmology contribute in an important way to the 
bispectrum. 

The second effect is a contribution to the bispectrum of 
local form which results from the curvaton mechanism. It 
is due to the transformation of entropy perturbations into 
curvature fluctuations which happens around the bounce 
point. Here, we focus on the second effect. 

A naive way to investigate the non-Gaussianity is to 
study the cosmological fluctuations to second order [7r| 
9 . Recalling the time evolution of the curvature pertur- 
bation given in Eq. (fTU|) and taking the next-to-leading 
order corrections into account, we are able to decompose 
the curvature perturbation into a Gaussian part and a 
non-linear one, where the second part evolves according 
to the following equation [83l - l86j : 



c 



NL 



W 2 )5s 2 - ^rV a Ss5s , (62) 



on long wavelengths. In the above formula, we have intro- 
duced the adabatic field a = J(<j) 2 + x 2 )^dt, the isocur- 
vaturc field s = (cf>x — X^)/ 17 , an d the trajectory angle 9 
satisfying cos = <p/& and sin# = x/a. 

As usual, we define a nonlinearity parameter Jnl to 
measure the amplitude of non-Gaussian perturbations. 
This parameter is defined by the relation 

3 



C = Cl + -JnlC 



L ■ 



(63) 



where Cl is the linear order value of curvature perturba- 
tion. As analyzed in the previous section, we consider 
the leading non-Gaussianity when modes are at super- 
Hubble scales, and thus the non-Gaussianity is of the 
local type[z3,|83|. 

We follow the time evolution of the curvature pertur- 
bations and use the solution (j2"Tj) . and then obtain the 
non-linear part of the curvature perturbation, 



NL 



H 



dt 



32m 



(1 



8d 2 Cm p i 



(64) 



where the main contribution is provided by the integral 
from t m to td during which the scalar x comes to dom- 
inate. Moreover, combined with the linear form of the 



Numerical analyses of inflationary preheating were performed 
with a massive scalar [68147CJ and a masslcss scalar l7ll [73l (see 
also [73|) respectively. 



9 The issue of inflationary non-Gaussianities has been studied in 
detail e.g. in Refs. [77t l78f | and comprehensively reviewed by 
[roll . Moreover, the non-Gaussian fluctuations during inflation- 
ary preheating have been analyzed in |67l ; see also [80-82] . 
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curvature perturbation (|56p. we obtain the leading term 
of the nonlinearity parameter 



5120 



N L 



d 2 C. 



(65) 



As discussed in the previous sections, there are two 
cases for the choice of d and C parameters depending on 
how the scalar field x becomes dominant. For Case 1, we 
get 

f NL ~ -3.3 , (66) 
by making use of C = (ir/Ad) 2 ; for Case 2, we obtain 



NL 



-5.3 



m 



d 2 M 4 



(67) 



In both two cases, we see that it is possible to generate 
a sizable non-Gaussianity in our model by this mecha- 
nism. Specifically, a large non-Gaussianity of local type 
can be obtained for a small value of d in Case 2. For 
example, taking the value of the parameters to be such 
that M = 100m which satisfies the theoretical bound 
discussed at the end of the Appendix, we find that with 
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-2 m 
M 



one gets Jnl — —5.3, a value which is con- 



sistent with current observational constraints which give 
—36 < Jnl < 100 at the 2a level, but a value which is 
not too small to be detectable in upcoming experiments. 

Although the mechanism discussed above is quite ef- 
ficient to generate a sizable non-Gaussianity, there are 
other sources which could potentially also give a large 
value of Jnl- For example, we expect that preheating 
effects will lead to a major contribution. In particular, 
they will affect the shape of the non-Gaussianities, a topic 
we did not study here. We will leave the detailed study 
of this issue to a future publication. 



bounce, a scenario in which the violation of the null en- 
ergy condition which is required to obtain a non-singular 
bounce is obtained via the addition of a ghost field. How- 
ever, the mechanism will also be effective in other real- 
izations of the matter bounce. 

As we have shown here, just before and during the 
bounce, the entropy fluctuations are amplified by an ef- 
fective tachyonic mass due to a kinetic resonance. We 
have called this process "kinetic amplification" . Whereas 
the scalar modes are amplified by this mechanism, the 
gravitational waves are not. Thus, our scenario leads to 
a suppression of the tensor to scalar ratio r. The kinetic 
amplification mechanism is more effective if the bounce 
is slower, since then the time interval during which the 
entropy fluctuations see an effective tachyonic mass is 
longer. The kinetic amplification mechanism is also more 
effective if the energy density in the entropy field initially 
dominates over the absolute value of the energy density 
in the ghost field. 

We also have taken a first look at non-Gaussianities 
from our mechanism. Specifically, we have discussed 
one source of local type non-Gaussianity generated in 
our curvaton mechanism. We find, for typical parame- 
ter values, the predicted amplitude is quite sizable to be 
observationally relevant. There are other sources, how- 
ever, which are expected to lead to even larger effects. 
Note that there have recently been some attempts to 
seed sizable non-Gaussianity in the framework of infla- 
tionary preheating [88l - [9lj . The distinguishing difference 
between their models and ours is that in our case the 
non-Gaussian fluctuations arise even for modes in the far 
infrared, modes relevant for current cosmological obser- 
vations. However, those in inflation are contributed by 
a narrow resonance which puts an upper bound on their 
wavelength (92l |. 
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VII. APPENDICES 

The the first part of this Appendix we study the evolu- 
tion equation for in the second we look at the evolution 
of the entropy fluctuations in the deflationary phase. 
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A. General Analysis of £ 



Consequently, its time derivative is given by 



As in the main text, we assume that matter is de- 
scribed by a collection of scalar fields. In this case, at 
linear order, the generic form of the equation for the time 
derivative of ( is 



H d>4 H 



(68) 



On large scales, the first term of the r.h.s of Eq. (f6"5)l 
is negligible. The second term describes the transfer of 
entropy to adiabatic fluctuations, the term we are inter- 
ested in. 

At first, let us see how the above multi-field formula 
applies to the double field case which is usually analyzed 
in the literature. For a model with two fields 0i and 02, 
we get 



H 0! <t>\ + 2 



2 <p\ 



(69) 



One notices that 



( yi v + ( ^ y = o 
K <t>\ + 4>1 4>\ + 4>l 



(70) 



and thus one obtains 

c = 



H 0, ^0? 



6o L 2 V W,? 



1 



^V 2 $ 
H 



2 l 0i 2 n 0? + 0| J ' 



(71) 



which is exactly what is obtained in the usual literature. 

In the following, we provide the detailed derivation of 
the equation (|68| for the time derivative of £ in the mul- 
tiple field model. We start with the standard perturbed 
Einstein equations, 

<h + AH$ + {2H + -m 2 )<S> = 4ttG6 P , (72) 
$ + H<$> = AirGSq , (73) 
-3iJ($ + #$) + V 2 $ = inGSp , (74) 



where 



sp = ]T[0\ 
<ty = Y^ 

i 

s P = y& 



0» + v ti 



(75) 
(76) 
(77) 



Recall that the gauge-invariant curvature perturbation 
variable ( can be expressed as 



2j 

C = $- -^($ + ii$) 

H 



(78) 



IT 

C = $ - ($ + - -^[4ttG(5p - H<$> 

TJ TT 

-3H($ + H<$>)] + -^-($ + H$) 

H HH 

= — r- UnGSp + 4ttG5o - V 2 $]H ^Sq 

H 2H 2 

= — ^[8ttG^ 0,(50; + 2F$-2V 2 $ + V 2 $] 

i 

HH 

+47rG-^- V 0^0i 

= H v 2 $ - 8ttG^- [V 0,(50, - V 0,(<50, + 3H6&) 
7? H . 

-Vy i 50 i ]+4 7 rG^V0 i 50 i 

XI 

= -|:V 2 $- 8ttG^ 53^*04 + 4^^^0,(50, 



H 

n 



V^ + SirGHY-^i- 



H 2H 2> 



^-^ghy 6 -^^)- 

H — 6: H 



(79) 



In the above calculation, we have used the background 
equation 



H = -iwGY' 



(80) 



a number of times. Moreover, it is convenient to intro- 
duce the Sasaki-Mukhanov variables as defined in the 
Section II, Qi = (50.; + Then we precisely obtain 

Eq. ©. 



B. Fluctuations in the Deflationary Phase 

Now we study the evolution of the entropy scalar x m 
the deflationary process. From studies of chaotic type 
inflation models, it is well known that the background 
field freezes out when the amplitude becomes of 
the order of the Planck mass, more specifically when 
Tn = 2 "' p ' . This leads to a nearly constant Hubble 

parameter H m = — ^ and an almost exponentially de- 
creasing scalar factor for the universe. This phase is cru- 
cial to achieving a nonsingular bounce, since the matter 
component, which effectively carries the negative energy 
density (namely, the Lee- Wick partner field or a ghost 
condensate), will grow rapidly and soon catch up to the 
density of 0. 

Assuming that the Lee- Wick partner -0 satisfies the 
following background equation of motion, 



i> + 3Hifj + M 2 tp = 



(81) 
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it is always oscillating around the stable point of its po- 
tential throughout the whole evolution of the universe. 
Before the deflationary phase, both the amplitudes ip and 
4> are growing in the same rate. Therefore, it is conve- 
nient to define a ratio 



where we have introduced an amplification factor C, 
which is determined by 



do = 



(82) 



K d 2 M> 



(86) 



which is almost constant until the moment t m , the end of 
the matter-dominated phase of contraction. Afterwards, 
solving the background equation of motion (|81[) , one ob- 
tains 



2m 



(t t m J 



(83) 



which evolves proportional to a -3 / 2 . One can compute 
that the contribution of ip catches up to that of 4> at the 
moment 



td t m ~t~ n 



7T m 

In 

2m d 2 M 



(84) 
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when 

We are able to compute the solution to the entropy 
field x as well. This yields x ^ a 3 - Then we continue 
to solve the perturbation equation in the infrared limit 
k — > 0, and obtain the amplitude of the isocurvature 
perturbation at the moment td as follows, 



- C 1 ' 2 5 Xm , 



4m 



{td — t m ) 



(85) 



indicating that the amplification on the isocurvature per- 
turbation in deflationary process is depressed by the ratio 
m/M but manifestly amplified by d 2 . The coefficient d 2 
is determined by the the ratio of two background scalar 
fields at low energy scale which is far away before the 
bouncing point. At that moment, the excitation of ip is 
strongly suppressed by the mass scale M as illustrated in 
the Lee- Wick model and the ghost condensate, and thus 
d 2 < m/M is satisfied commonly. 



Note that, the above analysis involves an important 
requirement that the coupling between the background 
scalar <p and the entropy field x is enough week, i.e., 
g < m/rripi. It is possible to choose a large value of the 
coupling constant which breaks this bound. This case 
corresponds to a stochastic resonance for the isocurvature 
perturbation before the bouncing phase, but would spoil 
the feature of scale-invariancc as observed by current ex- 
periments. A detailed analysis of stochastic resonance on 
entropy field around the bounce will be performed in an 
accompanied paper in a following up project. 
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